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Abstract 

The motility of spermatozoa of both Lytechinus pictus and Strongylocentrotus purpuratus sea urchin species is modulated by 
the egg-derived decapeptide speract via an oscillatory [Ca^^-dependent signaling pathway. Comprehension of this 
pathway is hence directly related to the understanding of regulated sperm swimming. Niflumic acid (NFA), a nonsteroidal 
anti-inflammatory drug alters several ion channels. Though unspecific, NFA profoundly affects how sea urchin sperm 
respond to speract, increasing the [Ca^^, oscillation period, amplitude, peak and average level values of the responses in 
immobilized and swimming cells. A previous logical network model we developed for the [Ca^^ dynamics of speract 
signaling cascade in sea urchin sperm allows integrated dissection of individual and multiple actions of NFA. Among the 
channels affected by NFA are: hyperpolarization-activated and cyclic nucleotide gated Na^ channels {HCN), [Ca^^- 
dependent CI" channels {CaCQ and [Ca^^-dependent channels {CaKQ, all present in the sea urchin genome. Here, using 
our model we investigated the effect of blocking in silico HCN and CaCC channels suggested by experiments. Regarding 
CaKC channels, arguments can be provided for either their blockage or activation by NFA. Our study yielded two scenarios 
compliant with experimental observations: i) under CaKC inhibition, this [Ca^^-dependent channel should be different 
from the Slol channel and ii) under activation of the CaKC channel, another [Ca^T channel not considered previously in the 
network is required, such as the pH-dependent CatSper channel. Additionally, our findings predict cause-effect relations 
resulting from a selective inhibition of those channels. Knowledge of these relations may be of consequence for a variety of 
electrophysiological studies and have an impact on drug related investigations. Our study contributes to a better grasp of 
the network dynamics and suggests further experimental work. 
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Introduction 

Fertilization is an important process in life. Reproductive 
success is attained by means of different strategies that increase the 
probability of gamete encounter. Several species, including sea 
urchins, produce spermatozoa with swimming patterns regulated 
by egg secretions. Strongylocentrotus purpuratus and Lytechinus 
pictus sea urchin spermatozoa swimming is modulated by speract, 
a decapeptide contained in the outer coating of the egg which 
diffuses in the sea [1,2]. When these sperm detect speract by 
means of receptors along the flagellum, an intracellular signaling 
pathway that regulates fluctuations of the intracellular Ca"'*" 
concentration ([Ca "^j) is triggered. The pathway involves 
production of cyclic nucleotides and alterations in the ionic 
permeability of the plasma membrane (Fig lA). These biochem- 
ical changes are timely translated to enhance the encounter of the 



sperm with the egg. Sea urchin spermatozoa swimming close to a 
surface describe circles, a very convenient trait to do imaging [3- 
5]. In the presence of speract gradients a cross-correlation of 
[Ca^'^ oscillations and path curvatures takes place [3,5-7]. 
Increases in the [Ca'^, are associated with sharp turning events 
(high path curvature) that are interspersed with periods of 
straighter swimming episodes (low path curvature). This swimming 
pattern is common to a wide variety of organisms with external 
fertilization [3,4,8-16]. 

Niflumic acid (NFA), a nonsteroidal anti-inflammatory drug, is 
able to block or modify several ion channels. Its lack of specificity, 
usually disadvantageous, turns out to be key to the profound effects 
it generates on how sea urchin sperm respond to speract. 
Immobihzed S. purpuratus sperm exposed to NFA respond to 
speract with [Ca^""^,; fluctuations that are larger, longer and with 
increased time intervals between them than control speract 
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Figure 1. Speract-activated [Ca^"*"] signaling pathway network model. A) Upper part: Schematic representation of the components of the 
signaling pathway triggered by speract in the sperm flagellum. Arrows traversing the membrane show ion fluxes. Arrows within the cell are indicative 
of causal relations. B) Bottom part: Signaling pathway operation diagram, black arrows correspond to activation, red lines to deactivation and yellow 
arrows can be activating or inhibitory depending on the relative state of the pathway elements being interconnected. Once speract binds to its 
receptor the several feedback loops are triggered according to the nature of the links involved. The concatenation of these loops leads to oscillatory 
stages of the whole pathway. The color code identifies corresponding upper and lower part components. Current models propose that the binding 
of speract to its receptor promotes the synthesis of cGMP that activates K+ selective and cyclic nucleotide-gated channels (KCNG) leading to 
membrane potential (V) hyperpolarization [3,4,7-11,18]. This V change first induces an intracellular pH increase via a Na + /H+ exchanger (NHE) 
activation, [18,51,52], stimulates hyperpolarization-activated and cyclic nucleotide-gated channels (HCN) [53-55], removes the inactivation of voltage- 
gated Ca^+ channels HVAand LVA [18,56] (CaV), and facilitates Ca-+ extrusion by Na + /Ca^+ exchangers (NCE) [51,52]. The opening of HCN and the 
influx of Na+ contribute to V depolarization, and concomitant increases in [Ca-'^]i and [A'a''"], further depolarize V. It has been proposed that the 
[Ca^^]- increases could lead to the opening of Ca- + -regulated Cl^ channels (CaCQ and/or Ca^ + -regulated K+ channels (CaKQ, which would then 
contribute to hyperpolarize the V again, removing inactivation from CaV channels and opening HCN channels, [3,4,1 8]. It is thought that this series of 
events is then cyclically repeated generating a sequence of V-dependent turns. B) Network model of the signaling pathway. The network can be 
envisaged as a circuit where each node represents an element of the pathway and links, either in the form of arrows or lines, correspond to 
connections determined in the bottom part of (A). The activating or inhibitory nature of the yellow lines depends on the value of voltage (V). Yellow 
nodes represent binary nodes (0,1), and the four brown nodes are ternary nodes that can take values 0, 1 and 2. Changes in the node states are 
determined by the connected nodes by means of a regulatory function (or truth table). As an illustration we present the case of the cGMP shown at 
the bottom left of (B). The first 3 columns in this table contain all the possible activation states of the cGMP regulators: GC, which is an activator; PDE, 
an inhibitor and cGMP (cGMP is a self-regulator); the fourth column shows the values for the function that correspond to each combination of the 
regulators. Additional nomenclature note: Speract receptor (SR); guanylate cyclase (GC); unknown Ca-^ channels sensitive to cAMP (cAMPCC); Ca^^ 
pump (CaP); dCa, dCI, dNa, dK are abbreviations for permeability changes in [Ca- + ], [Cl^], [Na + ] and [K+], respectively. 
doi:1 0.1 371/journal.pone.01 04451. gOOl 



responses [17] (Fig 2). Furthermore, in S. pwpuratus swimming 
sperm these alterations on flagellar [Ca^^, dynamics caused by 
NFA increase the speractinduced flagellar asymmetry resulting in 
more pronounced and prolonged sharp turns [4]. In L. pictus 
swimming sperm NFA inhibits chemotaxis (also in Arbacia 
Punctulata [15]), without altering their ability to detect the 
chemoatractant gradient [8]. 

It is because of the findings mentioned above that it is very 
important and interesting to unravel how NFA achieves these 
responses. There are not many tools that efficiently allow analysis 
of multiple variables regulating a transduction path; certainly a 
network model of a signaling pathway is one. In this regard using 
this approach to examine how NFA may achieve its remarliable 
effects on the speract signaling cascade seems warranted. This 
Boolean network model has already shown its ability to uncover 
and predict properties of this signaling cascade. In a previous 
publication we presented a theory-experiment investigation where 
the model predicted the possible participation of a CaKC channel 
and a voltage dependent Ca^^ channel. Experiments simulta- 
neously recording sperm [Ca ^, and swimming trajectories, and 



using specific blockers for these channels gave results consistent 
with the model predictions [18]. 

Analysis of multi-target drugs provides a means for revealing 
interrelations and synergistic outcomes among various compo- 
nents of a biochemical pathway, evidencing systemic behaviors. 
The complexity of the experimental realization of such a research 
is evident. Taking into account our previous work with NFA 
[4,8,17] we focus on this inhibitor as an example. In particular, 
there is experimental evidence that NFA inhibits [C a""*] -depen- 
dent [CI ] Channels (CaCC) and Hyperpolarization-activated and 
Cyclic Nucleotide-Gated Channels {HCN), and stimulates [Ca^""^- 
dependent [K"^] Channels (CaKC) of large conductance (BKca) 
[19-35]. However, it is not known it CaKC of small conductance 
(SKca) are blocked or activated by NFA [36]. In addition, 
blockage remains a possibility for other types of K""" channels such 
as K^j-p [37]. Hence, here we explore the consequences of 
assuming either activation or inhibition of CaKC channels. 

Here we use our model to explore the effect of NFA on the 
speract signaling cascade of S. purpuratus and L. pictus sperm, 
always keeping in mind previous experimental results. We focus on 
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Figure 2. Niflumic acid increases tKie [Ca^'*'], mean, maximum peai<, amplitude and interval between successive [Ca^'*'] fluctuations. 

Experimental determination of the Co- + -sensitive colorant Fluo-4 fluorescence in the flagellum, averaged along its length, of sperm exposed to 
100 nM caged speract, in the absence of NFA (black trace) and the presence of 10 /(M of NFA (red trace). Concentration is measured in florescence 
intensity (F.I) units and time in seconds according to the scale shown in the inset on the upper right part of the figure. The figure is a modified version 
of the equivalent shown in [4]. 
doi:1 0.1 371 /journal.pone.01 04451 .g002 



four individual-cell quantities measured previously by Wood et al., 
in [4], namely, average [Ca^^],- level, amplitude, peak and 
frequency. We use them to compare the model predictions 
between oscillations in the wild type (WT) untreated network and 
in the NFA treated case, considering the two scenarios, depending 
on wheter NFA inhibits or activates CaKC channels. 

Methods and Models 

The logical network 

The logical signaling network corresponding to the SASP, first 
introduced in [18], consists of nodes interlinked according to 
Fig lA, representing the principal components involved in the 
signaling cascade: ion channel activities, intracellular ion and 
molecular concentrations and the membrane potential amongst 
others. The network is shown in Fig IB and the nomenclature is 
explained in its figure caption. To analyze the dynamics of the 
network, we implemented a discrete formulation that is a 
generalization of the Boolean approach and that has proven to 
be revealing for the gene regulation dynamics of many systems 
[38-42], as well as other ceU signaling networks [43]. In this 
approach, the dynamical state of the network consists of a set of 
discrete variables {(Ti, 02, . . . , a,,}, each representing the state of 
a node. For this particular network, most of the variables take on 
two values, 0 and 1, depending on whether the corresponding 
element is absent or present, closed or open, inactive or active, etc. 
However, an accurate description of the dynamical processes in 
the network required four nodes to be represented by three-state 



variables: the membrane potential (hyperpolarized 0, resting 1, 
and depolarized 2); the low and high threshold voltage-gated 
Ca^^ channels (inactive 0, closed 1, and open 2); and the 
intracellular [Ca^^] concentration [Cfl^+], (basal 0, tonic 1 and 
supratonic 2). The state of each node a„ is determined by its set of 
regulators (which are some other nodes that also belong to the 
network). Let us denote as ff^j , (Jfii^ ■ • • 5 ^tiii the Ic regulators of (J/j, 
Then, at each time step the value of (T„ is given by 



a„it+\) 



(0, 



'"2 



it), 



(0 



(1) 



where F„ is a regulatory function constructed by taking into 
account the activating/inhibiting nature of the regulators. Each 
node has its own regulatory function. For the construction of these 
regulatory functions, we have made use of extensive biological 
knowledge, mainly of an electrophysiological nature, available to 
us in the literature and in our own laboratory. An illustration of 
such regulatory functions is given in Fig 1 for the cGMP node. For 
the case of [Ca"^], which is the main concern in our study, we 
have that it is regulated by 6 nodes (see origin of incoming arrows 
and red lines oidCa in Fig IB). LVA, HVA and cAMP-dependent 
Ca^**" channels are activators, i. e., their activation at time t, wiU 
favor and increase in [Ca^'*] at the subsequent time t + 1. There 
are also two inhibitory nodes: the Ca^"** pump and the Na''"/Ca^''' 
exchanger which when activated at time t will favor a decrease in 
[Ca^^] at time t + 1 . Finally the sixth node is the [Ca^'*] node itself 
that acts as an inactivator, hence favoring a decrease in its value at 
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the following time step. The balance among all the above inputs 
determined by physiological considerations is expressed in the 
form of a truth table with 432 rows, which constitutes its 

regulator)' function. In the supplementar)' material SI the 
regulatory functions for all the networks nodes are .shown 
explicitly. With this model we can therefore observe in silico the 
effect of altering certain elements relevant to the pathway. In this 
paper we consider the case of NFA-sensitive channels: HCN, 
CaKC or CaCC. In order to test the effect of NFA on the network 
evolution, we proceeded with the deletion or activation of these 
channels according to the two scenarios mentioned above, 
individually, by pairs {HCN and CaCC, HCN and CaKC, CaCC 
and CaKC) and all three simultaneously. Within our model, node 
elimination takes place by assigning a zero value to all the 
regulatory entries in which it participates during the whole 
network evolution. Hence, each channel wiU be closed, even if its 
regulators are activated, thus simulating the blocking effect of 
NFA. Activation is put into practice by keeping the affected 
channel permanently open with a value of one. We calculated the 
modeled speract-triggered [Ca^+], fluctuations averaging the 
[Ca^"*"] node value over 10' different initial conditions of the 
entire network (with all nodes present) and compared the result 
with equivalent calculations done on the treated networks. We 
repeatedly checked that after 10'^ initial conditions the averaged 
[Ca^+] values essentially coincided with the complete initial 
condition calculations. With the averaged intracellular [Ca^+] ion 
concentration, determined from the model dynamics, which we 
shall denote by [Ca^ + ],, we have a representation with a finer 
resolution, more adequate for comparison with continuous 
experimental measurements. 

Steady state analysis 

For this type of network (finite number of nodes which take a 
finite number of values), all initial conditions lead to a periodic 
behavior where the network configuration is rephcated after a 
certain number of steps. The time required to reach this condition 
is known as the transient time; it is important to mention that, 
independentiy of the initial condition, it is shorter than 45 steps. 
The attractor is reached after 22 iterations and the number of 
iterations between the repeated configurations is the period. These 
periodic solutions are the attractors of the network dynamics. For 
the wild type (WT) network (with all nodes present), 88.9% of all 
speract activated initial conditions lead to a period-4 attractor, 
while the remaining 11.1% converge to a period-8 attractor [18]. 
The total amount of initial conditions which reach a particular 
attractor is called the basin of attraction. In order to understand 
the differences between the averaged [Ca^"*^ concentration 
overline[Ca^^]^ steady-state behavior of NFA- treated network 
and the WT, we calculate modifications of the network attractors 
after alteration of all combinations of NFA-sensitive channels, 
paying attention to the ensuing number of attractors and their 
associated periodicit)'. We make the comparisons between the 
characterization of the periodic behaviors of [Ca^+],- and the 
period of the network attractors. This is done by means of a 
Fourier spectrum analysis of [Ca^+], from which we can 
determine differences in the temporal behavior of WT and 
NFA- treated network dynamics. We should point out that we have 
explored the effect of small perturbations on the regulatory 
functions, such as modifications in the outcome of any row in the 
truth tables of all the nodes. This holds particularly well for 
[Ca2+] and is a measure of the high robustness of the dynamics we 
are implementing and hence of the reliability of our results. 



Results 

Individual Channel Results 

In the first column of Fig 3, we show the [Cfl^ + ], fluctuations 
after a transient of 150 time steps, determined by the model as a 
function of time, for the original network and for the network with 
the channel specified in the figure blocked (A, D and G) and 
activated (J). We should point out that an extensive statistical study 
shows that steady state conditions are attained with less than 150 
iterations. Note the changes in the peak, mean and amplitude 
values, comparing the colored lines against the WT (black) curve. 
For a more quantitative comparison, in the second column (B, E, 
H and K), we show the bar plot of the average over 50 time steps 
of the amplitude and mean of the series shown to the left. In the 
third column (C, F, I and L) we focus on the temporal 
characterization of the series by means of the Fourier power 

spectrum of the [Ca^+J,- fluctuations. For each row, the channel 
altered is the same. Looking at the right column, the lines indicate 
the amplitude of the dominant frequencies in th(^ Fourier 
decomposition of the [Ca^+]; fluctuations depicted on the left 
column. The magnitude of each line indicates the weight of an 
underlying periodic contribution to the [Ca^+], time series of a 
given Fourier decomposition mode. For the WT case, in Fig 3C, 
the first black line to the left, at the frequency 1/8 is the 
fundamental period-8 Fourier mode and reflects the contribution 
to the [Cfl^+], of the period-8 attractor obtained from the network 
dynamics. The third line at 3/8 is its harmonic. The middle line at 
frequency 1/4 is the result of the averaging procedure of the 
[Ca^"*"] node network dynamics pc'riod-4 attractor together wtih 
the second harmonic of the period-8 attractor. We shall say that 
the [Ca^+], time series has a richer temporal behavior if its Fourier 
decomposition is more complex. Blockage of the HCN channel in 
our network model produces a more elaborate temporal behavior, 
which is registered in the Fourier power spectrum by the 
appearance of multiple lines that correspond to different modes 
in the Fourier decomposition. An opposite effect is observed with 
the blockage of the CaCC channel, where the Fourier decompo- 
sition consists of only one red line at the frerjuency value of 1/4 
that corresponds to a unique period-4 attractor of the [Ca^"*"] 
dynamics. For the case of CaKC blockage the Fourier spectrum 
remains unchanged with respect to the WT; this is consistent with 
the persistence in the blocked network of the WT attractors. 
Finally, activation of the CaKC channel produces 3 peaks at 1/7, 
2/7 and 3/7 frequency values which are related to a period-7 
attractor obtained by the network dynamics. 

Double Channel Results 

With regard to the temporal behavior, a similar analysis to the 
one performed for Fig 3, reaffirms the trends found from the 
individual channel alterations (Fig 4), namely: i) HCN blockage 
leads to a more complex Fourier Decomposition, consistent with a 
more elaborate behavior of the [Cfl^ + ] ■ series shown on the left, ii) 
CaCC blockage reduces the Fourier components leading to a 
simpler periodic behavior, iii) CaKC activation enhances the effect 
of the concomitant deletion of either HCN or CaCC. 

Triple Channel Results 

With the joint alteration of the HCN-CaCC-CaKC channels the 
differences of the effect of blocking or activating the CaKC channel 
can be appreciated in Figs. 5 and 6 respectively. Fig 5A is a 
segment of the [Ca^ + ]; time series once steady-state conditions 
have been reached. Fig 5B is the Fourier spectrum calculated over 
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Figure 3. [Ca^'*'] dynamics comparison between WT networl< and blocl<ing an individual NFA-sensitive channel. The [Ca-+] steady 
state fluctuations calculated form averaging over 100,000 initial conditions. For each initial condition, [Ca-+],. will take value 0 if it is in the basal state, 
1 if it corresponds to a tonic state and 2 if it is supratonic. After averaging it will take values within the range [0-2] with a resolution of 10^. 
Concentration units are arbitrary, they comply with the above restrictions and are set for comparative purposes, time is measured by simulation steps 
and frequency refers to the fraction of a cycle covered by a simulation step, i.e cycles per simulation step. The values of [Ca-+]^ determined by the 
logical signaling network model, in the above mentioned A.U. (Arbitrary Units), are shown in black for the wild type case and in colors for the network 
with deletion of the node indicated in the figure: In the first column note A) for HCN, the loss of regularity; D) For CaCC, the increase in amplitude, and 
G) for CaKC on scenario 1, a higher average and peak values. J) CaKC on scenario 2 (activating CaKQ, the decrease in peak, average and amplitude 
values. The central column (B, E, H and K) is the comparison between amplitude and mean of the [Cfl^+], dynamics between the WT curve and those 
with altered channels. The colors are the same as in the first column. Right column (C, F, I and L): Fourier Power spectra obtained by the curves on the 
left. 

doi:1 0.1 371 /journal.pone.01 04451 .g003 



1000 steady-state time steps. Besides tlie coincidence of the WT 
period 4 and 8 peaks with their harmonics as before, a new feature 
that arises in the dynamics of the blocked system is the occurrence 
of a period-9 peak and one of its harmonics. This is an indication 
of a richer temporal behavior. If we perform a running average of 

the [Ca^ + ], series over a window of 4 steps, the smaller 
fluctuations are dampened and the behavior of an envelope at 
larger scales becomes evident. This is shown in Fig 5C where a 
recurrent module of 72 time steps comes to light. Note that 72 is 
the minimum common multiple of the period 8 and 9 Fourier 
components, this is a manifestation of a period locking phenomena 

at the level of the [Ca^ + ], network attractors. If the running 



average is performed now over an 8 step window, then the WT 
[Ca^+], fluctuations are completely leveled out, as weU as the 
period-8 component of the blocked case network, allowing for the 
period-9 component to surface (see Fig 5D). When activation of 
CaKC channel is considered, we encounter the behaviors shown in 
Fig 6 which are determined as in Fig 5. In Fig 6A the emergence 
of a 3 step module absent in the WT is noticeable. This new 
feature is reflected in the period 3 Fourier component shown in 
Fig 6B. When a 4 step average is performed, the WT is strongly 
smoothened and a period 8 structure emerges while the structure 
of NFA treated case remains similar (see Fig 6C). In Fig 6D we 
show that though the WT is completely flatten to its over-all 
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Figure 4. Effect on the Cal cium d ynamics of altering the different NFA-sensitive channels taken by pairs. Average over 100,000 initial 
conditions of the steady-state [Cfl- + ], dynamics for the WT netw/ork is shown in black as in Figure 2. The columns are also distributed as in figure 2. 
A), B) and C) Blocking HCN and CaCC channels case (green line); D), E) and F) Blocking HCN and CaKC (dark green); G), H) and I) Blocking CaCC-CaKC 
(turquoise); J), K) and L) HCN blocked and CaKC activated; M), N) and O) CaCC blocked and CaKC activated. Notice that cases with HCN blocked 
produce non-regular [Ca^ + ] dynamics, opposite to the regularity generated by the CaCC blockage. Elimination of CaKC (scenario 1 in text) generates 
an elevation in [Ca^+] concentration compared with blockage of CaCC or HCN. Activation of CaKC (scenario 2) generates a decrease in [Ccr-+], and 
peak, but the temporal behavior is more elaborate. 
doi:1 0.1 371 /journal.pone.01 04451 .g004 



average, under an 8 step running average, when the effect of NFA 
is taken into account a module 3 structure is preserved. For the 
converse case of averaging over 3 steps, it is the treated case that 
gets flatten while the 8 module consisting of two alternating 4 
modules structure of tlie WT survives. 

Overall Picture 

Given our in silica results, we can conclude from our data, the 
following: Blockage of the HCN channel affects mainly the 
temporal behavior of the [Ca^ + ] dynamics. This became evident 
in figures 3 and 4, where changes in the temporal structure of the 
Ca^+ curve as well as the number of Fourier components can be 
observed. Another visible effect is the reduction of the [Ca^ + ], 
amplitude in all cases where HCN participates. When CaCC is 
blocked, the main effects are the increase of the amplitude in the 



[Ca^ + ], curve and a recovery of the regularity. Fourier components 
are reduced in all cases in which CaCC is altered. Regarding 
scenario 1 , in which CaKC channel is inhibited, the principal effect 
is an elevation in the [Ca^+],- mean and maximum peak. In scenario 
2, after activating the CaKC channel, the overall effect is quite 
similar to the one observed by HCN elimination: reduction in 

[Ca^ + ], amplitude, peak and mean but, an increase in Fourier 
modes. Though all of these effects are reproduced under 
combinations of the above channel alterations, it is possible to 
establish the dominance of one channel alteration over another. For 
example, with regard to the fluctuation amplitude we have that the 
effect of deleting HCN overrides CaCC blockage, however, there is 
a recovery of regularity, situation in which the effect on temporal 
behavior of CaCC blockage overrides HCN elimination (Fig 4A). 
The same temporal observation is obtained by looking at the 
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Figure 5. Effect of deletion of tKie three NFA-sensitive channels corresponding to scenario 1. A) Steady state time evolution of [Ca^+], 
time series calculated as described in Fig 2. Average oscillations for the wild type are in black while oscillations under deletion of the CaCC-CaKC- HCN 
nodes are in red (the NFA-blocked network). For the wild type, two contiguous 4 element modules can be identified that together constitute a 
recurrent 8-step module. For the treated network different 4-step and 8-step modules are indicated in the figure. B) Fourier spectra calculated from 
1000 steady-state steps of the [Co-+] ■ time series shown in (A). Period 4 and 8 Fourier modes and their harmonics are shown in black. For the NFA- 
blocked network the spectrum shown in red, determined from 1000 points of the steady-state [Ca- + ] time series, is richer due to the appearance of a 
period 9 Fourier mode with its harmonics, besides the period 8 contribution. C) A four element running average of the series corresponding to (A) 
including the initial transient time manifests an underlying repeated 8 module for the wild type and 72 module for the NFA-blocked network case. 
This last module is the minimum common multiple (MCM) of the period 8 and 9 Fourier modes shown in (B). For the wild type case, a period-8 
module surfaces which is the MCM of the period 4 and 8 Fourier modes shown in (B). D) Result of performing a running average as in C) using an 8 
element window instead. Here the wild type oscillations are wiped out (black graph) while a 9 module envelop (in red) is evidenced for the triply 
blocked NFA treated case. The insert is an amplification that shows this behavior in detail. 
doi:1 0.1 371 /journal.pone.01 04451 .g005 



Fourier power spectrum for the case of the joint blockage of HCN 
and activation oiCaKC. Both alterations result in an increase of the 
Fourier modes (Fig 4K) not only with regard to the WT but also to 
their action taken separately (Figs. 3C and 3L). This last result has a 
physiological explanation, because the effect of blockage of a 
channel which allows a cation entrance {HCN) or the activation of a 
channel which allows the cation extrusion [CaKC) will have a 
similar result in the regulation of membrane potential: shorter 
depolarizations . 

Finally, the effect of altering the three channels together is 
analyzed. Here, we encounter the two scenarios mentioned 
previously, depending on whether CaKC channel is inhibited 
(Scenario 1, Fig 5) or activated (Scenario 2, Fig 6). 

Overall our main findings are the following: 

With respect to scenario 1, our analysis shows that only when 
NFA blocks simultaneously the three channels under consideration. 



we recover the experimental observations of [4] with regard to WT 
conditions for the model [Cfl^+], time steps, namely: 

a. The mean taken over the time increases (Fig 5A, C and D). 

b. The amplitude grows (Fig 5D). 

c. The peak values are higher (Figs 5 A, C and D). 

d. The time interval between oscillations increases (Fig 5B). 

In scenario 2 (Fig 6), the effect observed is a decrease in 
amplitude, period, and maximum peak, a trend opposite to the 
experiments. This behavior suggests the need of another [Ca^+] 
channel in the signaling pathway. An integrated presentation of 
our findings is condensed in Table 1 . The red arrows correspond 
to those cases in which CaKC channel is activated (scenario 2). 
Table 1 clearly shows that within our modeling, experimental 
results from [4] can only be recovered under the joint NFA 
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deletion of the three HCN, CaCC and CaKC channels. The first 
three columns of Table 1, state the cause-effect connections of 
individual channel deletions: the blockage of the CaCC channel 
increases [Ca^+],- amphtude fluctuation, the blockage of the CaKC 
increases the [Ca^+],- peak and mean values, and the blockage of 
the HCN channel confers a more elaborate temporal behavior. 



Discussion 

Summary 

We developed a novel analysis of the effect of the multi-target 
drug niflumic acid, based on a logical network model of the 
speract-activated signaling pathway. We selectively modified the 
functionality of three ionic channels which are sensitive to NFA: 



Table 1. Effect of altering the NFA-sensitive channels in the Mathematical Model. 
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correspond to scenario 1, in which the CaKC channel is inhibited. The bolded arrows indicate the alterations corresponding to scenario 2, with activation of CaKC 
channels. 
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Figure 7. Schematic representation of tlie effect of deletion or activation of CaKC cKianneis in the speract activated signaling 
pathway. A) the membrane potential dynamics. B), the [Ca^+] dynamics. For (A) the membrane potential is depicted in black for the WT network 
after speract activation. Hyperpolarization (lower than -40 mV) due to a [A^^], efflux via the KCNG channel, and the consequent opening of the 
voltage-dependent HCN channel which in turn depolarizes the sperm flagellum are shown. Depolarization (higher than -40 mV) opens LVA and HVA 
[Ca-+] channels. The increase in [Co- + ],. enhances depolarization. This last entrance of [Ca- + ] opens the CaCC and CaKC channels with the 
corresponding influx of [C/^], and efflux of [if that hyperpolarize the membrane potential. All these steps are cyclically repeated causing the 
[Ca^+] oscillation pattern depicted in (B). If the effect of the drug is an inhibition of CaKC channels, this would cause a bigger depolarization, because 
a decrease in the [if^], efflux. This is shown in the purple curve, notice that the amplitude is bigger than in the WT case. If the effect on CaKC is an 
activation instead, a bigger hyperpolarization is produced, due to the increase in the efflux of [K^].. Depolarization takes less time and is smalles than 
in the WT case for the same reason (blue curve). For (B), the WT [Ca- + ] dynamics is again depicted in black. Inhibition CaKC channels (purple curve), 
reduces the extrusion of [K^]j diminishing hyperpolarization hence delayning the closure of CaV channels. This will cause the elevation of [Ca^"*"], as 
well as the time between [Ca^+] peaks (the period). Overall, the activation of CaKC channels (blue curve), produces a shorter [Ca-+] oscillations due 
to the faster hyperpolarization, which in turn closes the CaV channels avoiding a big elevation of [Ca-'^]j. 
doi:1 0.1 371 /journal.pone.01 04451 .g007 



CaCC, HCN and CaKC. We calculated [Ca^+] mean, amplitude, 
maximum peak and periodicity and determined the conditions 
under which accordance with previous experimental findings are 
attained. Other approaches have been done in terms of alterations 
of NFA in sea urchin spermatozoa [8,44] . Although NFA is known 
to affect these three channels, we addressed the question of 
whether or not a concerted alteration takes place in the sea urchin, 
taking into account that most pharmacological results have been 
derived in mammalian channels [19,21,23,24]. It CaKC channel is 
inhibited, we concluded that the experimental observations of [4] 
are retrieved only under a concomitant action of NFA on all three 
channels. The degree of synergy of this action remains to be 
elucidated. If CaKC channel is activated by NFA, the need of an 
additional [Ca^+] channel comes to light. Furthermore, along our 
investigation, indications of the relative importance of individual 
channel deletions surfaced. Our model predictions call for the 
implementation of experimental protocols for their corroboration. 
Knowledge of this type of cause-effect relations may be helpful for 
the understanding of the actions of diverse molecular compounds 
and might contribute to elucidate operating mechanisms of drugs. 

Physiological Context 

The physiological effect that NFA may exert on the Ca^ + - 
activated K+ channel CaKC in the speract-activated signaling 
pathway (SASP), either activating or inhibiting, can be explained 
in both cases (Fig 7): 

i) When CaKC is blocked, changes in membrane potential are 
attained, the depolarization is more pronounced owing to the 
reduction of [K^] efldux (Fig 7A, violet line), hence, the 



voltage-dependent [Ca^+j channels (CaVs) will remain open 
for a longer period; consequently, CaKC blockage (scenario 1) 
produces a bigger [Ca^ + ] curve (Fig 7B, violet line). Since 
there is pharmacological evidence of the presence of Slo 1 -type 
Ca^ + -activated K+ channel in the SASP [18], and it is 
known that the Slol channel is activated by NFA [45,46], the 
above line of reasoning could hold if we consider the 
participation of a difiFerent CaCKC such as SKcn or another 
K+ channel in the pathway. Given the results of the last 
column of table 1, these findings support this alternative, 
hinted by the experimental discrepancy mentioned in the 
introduction. 

ii) In the case of scenario 2, after the first .R^CWG-dependent 
hyperpolarization, the activation of the CaKC will produce 
subsequent higher hyperpolarization values since this CaKC 
channel remains open. This in turn causes shorter lived 
depolarizations of a smaller magnitude, situation that leads to 
a premature closing of CaVs with the consequent decrease in 
[Ca^+],. Under the assumption that NFA activates the CaKC 
channel, given that NFA is a non-specific drug, the 
inconsistency of the results between the [Ca^+j, network 
dynamics with a CaKC channels activated in silica and 
experiments, favors the consideration of other [Ca^+] 
channels in the SASP sensitive to NFA, so far not present in 
our model. Suggestions along these lines are met by the 
sperm-specific CatSper channel, present in the sea urchin 
genome [47-49], which has recentiy been found to be a 
polymodal chemosensor in human sperm [50]. The perti- 
nence and validity of this last suggestion, is a subject for 
further study. 
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Our investigation is closely related to the temporal behavior of 
the [Ca^"*"]; fluctuations. A proper understanding of temporal 
properties of the dynamics is fundamental for the comprehension 
of the link between internal biochemical oscillations that regulate 
molecular motors in the flagellum and the oscillatory space 
exploration of gradients of sperm activating peptides present in sea 
urchin swimming. Timing issues related to this interchange appear 
to be of importance for fertilization. In particular our study of the 
effect of NFA alterations on time-related characteristics may 
provide important issues on this link. Experimental analysis 
motivated by this work on the effect of NFA in the chemotactic 
response triggered by speract is provided in [8] . Finally, we would 
like to point out that the approach we have developed here can be 
used as a platform for the study of other biochemical signaling 
regulatory networks. 

Supporting Information 

Material SI Regulatory tables for the Network Dynam- 
ics of the SASP. Each one of the 21 text files which compose this 
compressed file, contains the name of the element (the file name), 
the name of its regulators, the whole combination of the dynamical 
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